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ABSTRACT: Stringent limits on the interactions between dark matter (DM) and the stan¬ 
dard model (SM) can be set by studying how initial-state or final-state particles recoil 
against missing transverse energy (px)- I n this work, we improve, extend and correct LHC 
constraints on the interactions between DM and top quarks that are mediated by the ex¬ 
change of spin-0 s-channel resonances. A comparison of the LHC run-1 sensitivity of the two 
main search channels is presented, which shows that mono-jet searches are typically more 
restrictive than the -\-tt searches. We furthermore explore the reach of the 14 TeV LHC. 
The collider constraints are compared to the restrictions arising from direct and indirect de¬ 
tection as well as the DM relic abundance, and we also reflect on effective field theory (EFT) 
interpretations of the LHC exclusions. 
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1 Introduction 

Searches for signatures represent one of the main focus of the ATLAS and CMS collab¬ 
orations in their hunt for physics beyond the SM, because of their possible connection to 
DM. These searches can be classified based on the type of SM particles which recoils against 
the DM pair. In LHC run-1, ATLAS and CMS have examined a variety of signatures 
involving jets of hadrons, gauge bosons, top and bottom quarks as well as the Higgs boson 
in the final state (for a recent review of the experimental status see for instance [1]). 

The existing LHC studies are in most cases performed in the context of an EFT which 
describes the physics of heavy particles mediating the interactions between DM and SM 
fields, assuming that the mediators are heavy enough so that they can be removed as active 
degrees of freedom. It has however been realised early on [2, 3], that an EFT description 
of signatures is deemed to fail if the masses of the mediators are within kinematic 
reach, which can cause the sensitivity of the LHC searches to change significantly. In order 
to correctly account for both off-shell and on-shell effects in DM pair production different 
simplified models have been put forward, in which the contact interactions present in the 
EFT are resolved into single-particle s-channel or f-channel exchanges. By specifying the 
spin and gauge quantum numbers of DM and the mediators and requiring the interactions 
to be minimal flavour violating (MFV) [4], the parameter spaces remain low dimensional, 
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which in turn allows for a simple translation of bounds between experiments and theories 
(cf. [5] for a comprehensive overview on simplified DM models). 

In the present work we focus on the DM pair production from quark or gluon initial 
states, where the production proceeds via the exchange of spin-0 s-channel mediators. 
Under the assumption of MFV, the most relevant DM-SM couplings in this class of simplified 
models are those that involve top quarks. Two main strategies have been exploited at the 
LHC to search for scalar and pseudo-scalar interactions of this type. The first possibility 
consists in looking for a pp —> $ T + j signal [6], where the mediators that pair produce 
DM are radiated from top-quark loops, while the second possibility relies on detecting the 
top-quark decay products that arise from the tree-level reaction pp —^ $T + tt [7, 8]. Further 
progress in characterising the LHC signatures associated to DM top-quark interactions has 
been made in [9-12] for the mono-jet signal and in the articles [11, 13-16] for final states 
involving top-quark pairs. 

Our goal is to refine, to extend and to correct existing LHC constraints. As a first 
benchmark, we examine the constraints on the EFT that stem from LHC run-1 mono-jet 
searches as well as the different $T + tt channels with di-leptonic [14], single-leptonic [15, 16] 
and fully hadronic [16] top-quark decays. This exercise serves not only as an independent 
cross-check of the $T + tt analyses performed by ATLAS and CMS, but also as a validation 
of our Monte Carlo (MC) chain. We find that at present the strongest constraints that the 
LHC can place on effective DM top-quark interactions arise from mono-jet searches, and 
that this strategy is expected to remain the most powerful one also at future LHC runs. 
Our EFT results at both 8TeV and 14TeV are compared to the exact exclusions limits 
in the simplified DM models. This allows us to determine under which circumstances an 
EFT interpretation of the collider bounds is justified. By scanning the parameter space of 
the simplified models, we furthermore show that in both the scalar and the pseudo-scalar 
case the ATLAS and CMS searches cannot presently exclude parameters arising from purely 
weakly-coupled theories. As far as a comparison is possible, this Ending agrees qualitatively 
with the conclusions drawn in [11, 12], Our analysis reveals in addition that the + j 
searches generically exclude more parameter space than the + tt searches. We finally 
discuss the interplay of the various DM searches, including direct and indirect detection as 
well as the constraints from the observed relic abundance. 

The outline of this paper is as follows: in Section 2 we discuss the structure of the 
simplified DM models and the corresponding EFTs. This section contains in addition the 
formulas needed to calculate the DM-nucleon scattering cross sections as well as the DM 
relic density. The results of our phenomenological analyses of the DM top-quark interactions 
are presented in Section 3. In this section we discuss the present and possible future bounds 
that result from the different LHC search strategies, comparing the obtained limits to those 
arising from direct and indirect detection as well as the requirement not to overdose the 
Universe. In this context, special attention is payed to the differences in the results and the 
conclusions drawn when the calculations are performed in the simplified model framework 
or the EFT. We conclude and provide an outlook in Section 4. Additional material that 
might be of particular interest for the practitioner is relegated to Appendix A. 
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2 Preliminaries 


In the following we introduce the simplified models for the DM-SM interactions that re¬ 
sults from the exchange of a colourless scalar or pseudo-scalar mediator (Section 2.1) and 
discuss the relevant operators in the corresponding EFT (Section 2.2). To make our article 
self-contained we furthermore collect the formulas necessary to calculate the DM-nucleon 
scattering cross sections and the DM relic abundance (Section 2.3). 

2.1 Simplified models 

The relevant interactions between DM and quarks involving the exchange of a colourless 
scalar ( S ) or pseudo-scalar (P) mediator are parameterised as follows 

£ 15 5 dm (xx) S + gluYl ~ (w) S + * 5 dm (X75X) P + * 5 sm ““ (™) P ’ 

9 9 

where the sum is over all quarks and v — 246 GeV denotes the Higgs vacuum expectation 
value. In writing (2.1) we have assumed that the couplings of the mediators to quarks are 
proportional to the associated SM Yukawa couplings. This is motivated by the hypothesis 
of MFV, which curbs the size of dangerous flavour-changing neutral current processes [4] 
and automatically leads to a stable DM candidate [17]. Notice that the MFV hypothesis 
allows the mediator quark couplings to be scaled by separate factors g d ' and g u ’ for 
down-type quarks and up-type quarks, respectively. For simplicity, we have ignored this 
possibility when writing C and choose the same scaling factors g d ' P = gu’ P = g f° r a H 
quarks. While the DM particle x hi (2-1) is understood to be a Dirac fermion, extending 
our discussion to Majorana DM or the case of a complex/real scalar is straightforward [6]. 
In order to avoid the severe experimental bounds from the electric dipole moment of the 
neutron (cf. [18]), we take the spin-0 mediators S,P to be CP eigenstates and in addition 
assume that the couplings g and g^ are all real. 

Further constraints on our simplified DM models can in principle arise from existing 
and future LHC resonance searches in it final states. Including the one-loop process gg —> 
S,P —> it, one finds [10] that for weakly-coupled models the total it cross section is changed 
by only 0(1%). Such small effects are likely to remain unnoticed given that the theoretical 
uncertainty on the total it cross section is at the level of 5% at the LHC [19]. A di-jet signal 
arises in the simplified models (2.1) first at the two-loop level via gg —>• S, P —> gg. The 
strong loop suppression renders the contributions of S, P exchange to di-jet production 
unobservable at the LHC [10]. Since the SM portion of the Lagrangian (2.1) is not a 
electroweak singlet additional restrictions also stem from the fact that the mediators S, P 
necessarily have portal couplings involving the Higgs field. The resulting modifications in 
Higgs phenomenology are, however, model dependent and we do not study them in what 
follows. If the mediators have weak-scale masses, the couplings in (2.1) to light quarks are, 
to the best of our knowledge, unconstrained by direct and indirect collider searches that do 
not involve large amounts of For Ms,p < 10 GeV important constraints can however 
arise from quark flavour physics [20]. 

The signal strength in DM pair production does not only depend on the couplings 

S p s p 

5d’m anc ^ 5 sm an< ^ masses rn x and M$,p, but also on the total decay widths T$,p of the 
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mediators S, P. In the case of the scalar mediator, one finds the following results for the 
partial decay widths (see e.g. [ 10 ]) 


T{S^XX) = 




4ttj2 \ 3/2 

li) ms-2m x ), 


-> qq) — (si m ) 

F(S —>■ gg) = (g § M ) 


3 ml M s 
8irv 2 


27 v 3 v 2 Mg 





3/2 


9 (Ms - 2 m q ), 

4 m 2 \ 

M 2 ) ’ 


( 2 . 2 ) 


where 9(x) denotes the Heaviside step function defined by 9(x) = 0 for x < 0 and 9(x) = 1 
for x > 0 , while 

Fg(x) = 1 + (1 — x) arctan 2 ( ) . (2-3) 

Wx - 1 / 

The analogue expressions for the pseudo-scalar mediator are obtained from (2.2) by the 
replacements S — > P and 3/2 —>• 1/2 in the exponents, and the relevant form factor reads 

Fp(x) = arctan 2 ( ) . (2-4) 

Wx- 1/ 

At the loop level the mediators can decay not only to gluons but also to pairs of photons 
and other final states if these are kinematical accessible. The decay rates T(5 —> gg) and 
T(P -> gg) are however always larger than the other loop-induced partial widths, and in 
consequence the total decay widths T 5 and T p are well approximated by the corresponding 
sum of the individual partial decay widths involving DM, quark or gluon pairs. Notice 

g p s p 

finally that if Ms,p > 2 mt and g > #dm’ total width of S', P is dominated by the 
partial widths to top quarks due their large mass or Yukawa coupling. 


2.2 EFT description 

If the mediator masses Mg t p are large compared to the other scales involved in a given 
process, one can describe the underlying partonic reaction by means of an EFT. Integrating 
out the scalar and the pseudo-scalar mediator then gives rise to 

°s = jxxxqq, o g p = ™^xi5xqi5q, ( 2 . 5 ) 

at tree level as well as contact terms consisting of four DM or quark fields. In the case of 
the scalar operators O q s the suppression scale Ag is related to the mediator mass Mg and 
the fundamental couplings and g§ M by 


As 


vM 2 s \ 1/3 
9sm9dm/ 


( 2 . 6 ) 


and an analogous expression with S —> P holds for the pseudo-scalar operators Op 
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Figure 1. Examples of Feynman diagrams with an insertion of Og that contribute to a p T + j (top 
left), $T+tt (—» 2 jbb) (top right), fip+tt (—» jblv) (bottom left) or a $ T + tt (—>■ bl~vbl + v) (bottom 
right) signal. The operator insertions are indicated by yellow blobs, while the regular SM vertices 
are represented by black dots. 


Integrating out the top quark generates an effective interaction between DM and gluons. 
At the one-loop level, one obtains 


o G = 


CXc 


12vrA| 


XXG%G 0 ^, 


o~ = as 

G 8nAl 


X'TSXG^G^, 


(2.7) 


by employing the Shifman-Vainshtein-Zakharov relations [21]. Here G“„ denotes the gluon 
field strength tensor and G a ^ v = 1/2 e^ wXp G^ p its dual. At the bottom-quark threshold 
and the charm-quark threshold one has to integrate out the corresponding heavy quark by 
again applying (2.7). Note that this matching procedure is crucial to obtain the correct 
DM-nucleon scattering cross section associated with effective scalar DM-quark interactions. 

Constraints from mono-jet searches on the scalar and pseudo-scalar DM-quark interac¬ 
tions Og p involving the light quark flavours as well as the gluonic operators Oq have been 
discussed in detail in [9, 22 ], while analyses of the bounds on O b s P arising from + bb 
final states have been carried out in [7, 8 , 13, 16]. We instead focus on the effective in¬ 
teractions Og p containing top quarks. In Figure 1 we show representative graphs with 
an insertion of Og corresponding to the different dedicated search strategies that have so 
far been exploited to constrain DM top-quark interactions at the LHC. From top left to 
bottom right these are mono-jet searches [ 6 , 9, 11, 12], + tt with fully hadronic top- 

quark decays [16], + it where one top quark decays hadronically and the other one 
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semi-leptonically [ 8 , 11, 15, 16] and ftp + it with di-lepton final states [11, 14]. Future 
prospects and opportunities of the searches for DM heavy-quark interactions have been 
discussed in [10, 13]. 


2.3 DM-nucleon scattering and relic density 

In the case of the scalar operators 0| the cross section for elastic Dirac scattering on a 
nucleon is spin-independent (SI) and given by 



™Ld m NfN 

ttA 6 s 


( 2 . 8 ) 


where ?n re d = m x mp[ / {m x + mpr) denotes the reduced mass of the DM-nucleon system, 
m,]\r — 0.939 GeV is the average nucleon mass and /A ~ 0.30 (see [23] for a recent evaluation) 
is the effective DM-nucleon coupling. For the pseudo-scalar operators Op, on the other 
hand, the DM direct detection cross section is spin-dependent and momentum-suppressed 
by Existing direct detection experiments are hence not sensitive to effective pseudo¬ 

scalar DM-quark interactions. 

In the presence of the effective interactions (2.5) two different annihilation channels 
contribute to the total annihilation cross section. Tree-level annihilation into quarks will 
be dominant for m x > rrit, while annihilation into gluons via heavy-quark loops can give a 
relevant contribution for lower DM masses. Performing an expansion in the DM velocity v x , 
the total annihilation cross section for Og and Op take the form {av x ) s = b s v\ + 0{v§ 
and (av x ) p = ap + 0(v x ) where (see for instance [ 6 ]) 


b S = Yl 
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3rriqm x 

8 ttA| 
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^ 3m 2 m 2 

^ 2ttA 6 p 
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m q ) + 


a: 


27T 3 Ap 



(2.9) 


and the analytic expressions for the form factors F$(x) and Fp{x) have already been given 
in (2.3) and (2.4), respectively. Note that the total annihilation cross section associated 
to O g s is p -wave suppressed, while DM annihilation proceeds via s-wave for Op. 

Since 7 -rays are an unavoidable product of hadronisation, it follows from (2.9) that for 
pseudo-scalar interactions indirect detection experiments can provide relevant constraints 
on the parameter space. The corresponding velocity-averaged DM annihilation cross section 
is given by {&v x ) p — ap/2. Another indirect constraint results from the requirement not to 
overdose the Universe. In terms of the coefficients (2.9) and the observed DM abundance 
Q x h 2 ~ 0.11 [24], the predicted DM relic density can then be expressed as [25] 


, 2 x 3.2 • 10 -8 GeV” 2 2 fr > u2\ 3.8 ■ 10- 9 GeV" 2 

{p, x h ) s — n x h , (n x h ) p — 


b S 


ap 


Q x h 2 


( 2 . 10 ) 
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3 Phenomenology 


In this section we study in detail the phenomenology of the DM top-quarks interactions 
induced by spin-0 s-channel exchange. After describing the main features of our MC sim¬ 
ulations of the individual fix signals (Section 3.1), we turn to the EFT and discuss the 
present (Section 3.2) and the possible future bounds (Section 3.3) on the suppression scales 
A s t p that result from the different search channels. We then investigate the accuracy of the 
heavy top-quark approximation for the case of the mono-jet cross sections (Section 3.4). 
The current (Section 3.5) and projected (Section 3.6) limits on the parameters space of the 
simplified models are examined subsequently. In this context, we discuss the complemen¬ 
tarity and the interplay of the various DM search strategies, including the constraints from 
direct detection, DM-induced 7 -ray emission from dwarf spheroidal satellite galaxies of the 
Milky Way and the relic abundance, and assess the quality of the EFT interpretations of 
$x searches at the LHC. 

3.1 MC simulations 

Our predictions for the mono-jet cross section are obtained using the POWHEG BOX [26] and in¬ 
clude leading order (LO) fixed-order contributions, parton-shower effects and hadronisation 
corrections (LOPS). The needed partonic one-loop amplitudes are taken from the MCFM [27] 
implementation of the process pp — > H/A+ j -7 t + t~ + j . which is based on the analytical 
results of [28] for the scalar Higgs case ( H) and [29] for the pseudo-scalar Higgs case (A). 
Our new MC implementation has been validated by calculating the partonic mono-jet cross 
sections both in the context of simplified models and the EFT, finding perfect agreement 
with existing numerical results [ 6 , 9]. To determine the cross sections for the different $ T +tt 
signals, we have implemented the Lagrangian densities (2.1) in FeynRules 2 [30], generating 
a UFO output [31]. The actual event generation is performed at LO with MadGraph 5 [32], 
Our MC chain has again been successfully validated against the results of previous stud¬ 
ies [ 8 , 13, 15, 16]. Parton-shower (PS) effects and hadronisation corrections have in all 
cases been included by means of PYTHIA 6 [33] and jets reconstructed using the anti-fc* 
cluster algorithm [34] implemented in Fast Jet 3 [35]. A detector simulation has not been 
performed, since even without it we are able to reproduce the relevant fix search results of 
ATLAS and CMS within errors. 

In contrast to the recent theoretical analysis [11] we do not multiply our LOPS predic¬ 
tions by a K factor to mimic the impact of next-to-leading order (NLO) corrections. In the 
case of the mono-jet signal this is motivated by the observation that the infinite top-quark 
mass limit is a bad approximation if the px cut on the jet is large and/or DM is heavy [ 6 ]. 
However, only in the case of mt —> 00 are the NLO corrections to pp —> H/A + j produc¬ 
tion known (see e.g. [36, 37]), while an exact 0(c4) calculation with resolved top-quark 
loops is at the moment unavailable. Until such a computation is at hand, we believe it is 
more conservative not to take K ~ 1.6 from pp —> H/A + j production and to apply it 
in the mono-jet cross section computation. In the case of pp —> H/A + tt production, on 
the other hand, the exact top-quark mass dependence is known at O(a^) already for some 
time (cf. [38-40]), and NLO effects turn out to be small, leading to K ~ 1.2. Given that 
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the NLO effects are not flat over the entire phase space and that the experimental cuts 
imposed in pp —> H/A ± it and pp —?• $r ± it are not identical, we again prefer to be safe 
and not to include a K factor in our results for the $r ± it cross sections. 

The predictions for all signals are obtained using MSTW2008L0 parton distribution 
functions [41] and the corresponding reference value for the strong coupling constant. In 
the case of mono-jet production, we define p = pr = pp = £Hr/2 and evaluate this 
scale on an event-by-event basis. Here pr an d pf denote s the renormalisation scale and 
factorisation scale, respectively, and Hr = \Jm\ x A Vr j 1 + PT,ji ■ The invariant mass of 
the DM pair is denoted by m xx and PT,ji corresponds to the transverse momentum of the 
hardest jet. In the case of the fix + it processes, we have instead employed the dynamical 
scale p = pr = pf = £ (rrq ± m xx /2). In order to assess the theoretical uncertainties 
that plague the calculated cross sections, we study the scale ambiguities by varying the 
parameter £ in the standard range [1/2,2], Numerically, we find that the predictions for 
the mono-jet cross sections calculated in this way vary in the ballpark of ±40%, while in 
the case of the ]fc T ± it processes, slightly smaller variations of around ±35% are obtained. 

3.2 Status of EFT limits 

In the following we list the various cuts and the values of the fiducial cross section (ofid) 
of each individual $r channel. This information will then be used to set limits on the 
suppression scales A s,p that appear in the effective interactions (2.5) involving top quarks. 

Mono-jet channel 

In order to derive the most stringent constraints from existing $r + j searches, we employ 
the latest CMS results [42], which make use of 19.7 ft -1 of 8TeV data. The relevant 
selection cuts are 

Pt.ji > 110GeV, |r/j, < 2.4, Pr,j 2 > 30GeV, \r]j 2 \ < 4.5, A 4>jij 2 < 2.5, (3.1) 

where A (f>j 1 j 2 is the azimuthal separation of the two leading jets, which are reconstructed 
using a radius parameter of R = 0.5. Another important selection criterion is the imposed 
jet veto [22], which rejects events if they contain a tertiary jet with PT,j 3 > 30 GeV and 
177y 3 1 < 4.5. The CMS measurement is performed in seven distinct $r regions, and we find 
that in the case of the operators O l s P the highest sensitivity is obtained for $r > 450 GeV. 
The corresponding 95% confidence level (CL) limit on the fiducial cross section reads 

Cfid(pp-t $ T + j) < 7.8 fb. (3.2) 

Fully hadronic $ T ± it channel 

The recent ATLAS search [16] looks for a $R-{-it (—>• 2 jbb) signal. In this analysis based on 
20.3 fb -1 of 8TeV data, jets are clustered with R = 0.4. To pass the trigger either five jets 
with pr,j > 55 GeV or four jets with px.j > 45 GeV one of which is identified as a bottom- 
quark ( b) jet are required. Events are only selected if they have at least five reconstructed 
jets, out of which two or more are 6-tagged, and they fulfil 

> 200GeV, \rjj\ < 2.5 , A^ > 1.6. 
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(3.3) 



Here b\ denotes the 6-jet with the highest transverse momentum. The 6-jet tagging efficiency 
is taken to be 70% here and in what follows. Based on this selection requirements, the 
ATLAS collaboration is able to set the following 95% CL limit 

<Aid (pp —> $t + 2j'66) < 2.0 fb. (3.4) 

Single-lepton $ T + it channel 

For what concerns the $p-\-it (—> jblu) mode, we again rely on the ATLAS results [16] (see 
also [43]), since this search turns out to be slightly more constraining than the dedicated 
CMS analysis [15]. As we have explicitly verified, comparable limits on the suppression 
scales Acan also be obtained by recasting the searches [44, 45] for top-squark pair 
production in the single-lepton final state. The trigger employed in [16] requires exactly 
one lepton (l = e, p) with pp,i > 25GeV and \r)i\ < 2.5 as well as four or more jets, where 
one jet is 6-tagged and all satisfy \iy \ < 2.5. Events are selected, if they pass the cuts 

> 270 GeV, pt j 1 > 80GeV, Pt j 2 > 70GeV, Pt j 3 > 50GeV, 

(3.5) 

PT,j 4 > 25 GeV, PT,b 4 > 60 GeV, m jjj < 360 GeV, A 4>f£ T > 0.6. 

Here mjjj is the three-jet invariant mass [43] and / = l. ji, j r 2 - Furthermore, the angu¬ 
lar separation between the lepton and the leading jet (6-jet) has to satisfy A R[j 1 < 2.75 
(A Ri bl < 3.0), the transverse mass my formed by ppi and $ T has to exceed 130 GeV and 

$ T j 'Y^ l=l PT,j n > 9 V GeV is required. The kinematic invariant arriT 2 [46-48] has to fulfil 
am,T 2 > 190 GeV. These requirements lead to the following 95% CL limit on the fiducial 
cross section 

&fid{pp —■► $t + jbl) < 0.5 fb . (3.6) 


Di-lepton $ T + tt channel 

We finally consider the results of the CMS search for a + tt (—> bl~Pbl + i') signal [14], 

performed on a 8 TeV data sample that corresponds to an integrated luminosity of 19.7 fb -1 . 
The basic selection requirements are pxj > 30 GeV, \rjj\ < 5, px,i > 20GeV, |r//| < 2.4, 
mu > 20 GeV, \rriu — 91GeV| > 15 GeV and the jet radius is R = 0.5. In addition, the 
following four cuts 

2 2 

$ T > 320 GeV , PT, 3 n < 400 GeV , ^ p T ,i n > 120 GeV , A0 Wa < 2 , (3.7) 

n= 1 72—1 

are applied to separate signal from background. The relevant 95% CL limit on the fiducial 
cross section is 

a M(pp * $t + bl~bl + ) < 0.15 fb . (3-8) 

Comparison of current constraints 

The two panels in Figure 2 show the 95% CL bounds on the suppression scales A 5 (left) 
and Ap (right) that derive from the individual search strategies discussed before. The 


- 9 - 




250 


250 


200 - 



5 10 


50 100 
i* [GeV] 


500 


200 - 



5 10 


50 100 
i* [GeV] 


500 


Figure 2. Lower 95% CL limits on the suppression scales A s (left) and A p (right) that derive from 
the mono-jet (red), $ T +tt (-* 2 jbb) (blue), $ T +tt (—» jblv) (green) and $ T +tt (—» bl~ubl + v) (yel¬ 
low) searches after LHC run-1. The widths of the bands reflect scale uncertainties. 


widths of the coloured bands illustrate the impact of scale variations. The corresponding 
relative uncertainties amount to around ±8% and ±6% in the case of the mono-jet and 
the ftp tt signals, respectively. We see that for both the scalar and the pseudo-scalar 
operator Og and Op, the bound (3.2) on the mono-jet cross section provides to the best 
constraints at the moment. Numerically, we obtain Ag > 145 GeV (A p > 160 GeV) for 
m x < 100 GeV, if theoretical uncertainties are included. We add that the LOPS mono¬ 
jet cross sections, on which our limits are based, are by roughly 45% smaller than the 
corresponding LO fixed-order results. A effect of similar size has been observed in the case 
of the operators Oq [22], which implies that the impact of the jet veto can be well modelled 
by working in the EFT (2.7), where both the mediators S, P and the top quark have been 
integrated out. 

Turning to the bounds arising from the $p + it channels, we first observe that our 
limits are in full agreement with the ones reported in [14, 16]. This shows indirectly that 
the reconstruction efficiencies are close to 100% in these analyses. It furthermore suggests 
that the signal reconstruction is independent of m x and does not depend on whether one 
considers the insertion of Og or Op. Second one sees that the strongest constraints stem 
from the single-lepton limit (3.6). This search allows to set a lower bound of A s,p > 120 GeV 
if DM is lighter than about 100 GeV. The corresponding limit for both the di-lepton and 
the fully hadronic pp+it channel amounts to A s,p > 90 GeV. Notice that in the case of the 
ftp + it modes the constraints on the suppression scale entering the operators Og and Op 
are very similar for light DM, while in the mono-jet case the limits on Ap are by roughly 
10% stronger than those on A g. This feature can be understood by observing that the 
ratio of the multiplicative factors appearing in the operators Oq and Oq reads 3/2 Ap/Aj|. 
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Figure 3. Projected lower 95% CL limits on A 5 (left) and A p (right) from mono-jet (red) and 
$ T + it (—> jblv) (green) searches at the 14TeV LHC. The shown predictions assume an integrated 
luminosity of 25fb~ 1 . Scale uncertainties are indicated by the widths of the coloured bands. 


From (2.7) one would hence expect that for a given mono-jet cross section the restrictions on 
A p are by a factor of (3/2 ) 1 / 3 ~ 1.14 better than the limits on A 5, and this is to very good 
approximation what one finds. Finally, realise that for m x > 100 GeV the limits on A 5 all 
fall off faster than the bounds on Ap. This property is related to the fact that the S —>• XX 
squared amplitude is proportional to — 4 m x , while for P —> XX one instead has 

3.3 Prospects of EFT limits 

It is also worthwhile to investigate how the reach on the suppression scales A g p might 
improve at 14 TeV. As we have seen the mono-jet channel and the + tt single-lepton 
mode provide at present the two most stringent constraints, and this situation is unlikely 
to change at future LHC runs. We therefore focus on these two search strategies, discussing 
the relevant experimental cuts and the associated SM backgrounds for each signal in turn. 

Mono-jet channel 

In the case of the mono-jet signal, we apply the event selection criteria that have been used 
in the sensitivity study by ATLAS [49]. They are given by 

PT,j 1 > 300GeV, | rj h \ < 2.0, p T j 2 > 50GeV, \q n \ < 3.6, A 4> jtr > 0.5, (3.9) 

and jets are defined with R = 0.4. Events with a third jet of pr,j 3 > 50 GeV and | rj J3 < 3.6 
are vetoed and > 800 GeV is employed, since out of the three thresholds discussed 
in the ATLAS study this cut provides the strongest restrictions. Notice that compared 
to (3.1) the Ptj'd PT,j 2 an d requirements are increased to avoid pile-up and to enhance 
the signal-over-background ratio. In order to determine the limits on A gp, we follow [49] 
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and take 


(3.10) 


°fi “{PP &T + j) = 5-5 flb, 

assuming a total systematic uncertainty on the SM background of 5%. 

Single-lepton $ T + it channel 

Our forecast for the $ T + it (—> jblv ) channel is based on the cuts introduced in the 
ATLAS benchmark study of top-squark pair production [50]. Specifically, we require a 
single lepton with pxi > 25 GeV, \r)i\ < 2.5 and in addition four or more jets with one b -tag 
all satisfying |rp \ < 2.5. Jets are found using R = 0.4 and we impose 

$x > 550GeV, PTji > 80GeV, Pt j 2 > 60GeV, Pxj 3 > 40GeV, 

(3.11) 

Px,j 4 > 25 GeV , 130 GeV < nrijjj < 205 GeV , A 4>f$ T > 0-8 , 

where / = j \, j?- The requirement on the transverse mass calculated from px.i and $x is 

m-T > 350 GeV and we ask for $x/\jY^n =l Pr,j n > 15 \/GeV. Compared to (3.5) the $ T 
and mx selections are significantly stronger in (3.11), which allows to better disentangle 
signal from background. For the above cuts, the total SM background amounts to [50] 

crf£(pp$ T + jbl) = 0.13 fb, (3.12) 

and has a total uncertainty of 7%. Notice that the cuts in (3.11) are not fully optimised 
for our purposes. The future constraints that we derive using them should therefore be 
conservative. 

Comparison of future constraints 

In Figure 3 we present our projection of the 95% CL limits on the suppression scales A s (left 
panel) and A p (right panel). In the mono-jet case, we observe that with 25 fb ^ 1 of data, cor¬ 
responding to the first year of running after the LHC upgrade to 14 TeV, one may be able to 
set a bound of A s > 205 GeV (Ap > 220 GeV) for m x < 100 GeV. Compared to the present 
mono-jet limits, this corresponds to improvements by a factor of 1.4. With 300 TV 1 of ac¬ 
cumulated data, we arrive instead at A 5 > 230 GeV (A p > 250 GeV). Collecting 3000fb -1 
will not allow to notably improve these limits, which shows that at 14 TeV the reach of 
the $x + j channel is not statistically limited, but limited by the systematic uncertainties 
associated to the imperfect understanding of irreducible SM backgrounds. Finding ways to 
overcome these limitations will be crucial to exploit the full physics potential of mono-jet 
searches to be carried out at later stages of the LHC. In the case of the fix + it search in 
the single-lepton final state, one observes that in the first year of data taking at 14 TeV the 
present bounds can be improved by a factor of 1.2 only. Such an improvement will allow to 
exclude scales A s,p > 140 GeV for m x < 100 GeV. Since the fiducial cross section (3.12) is 
compared to (3.10) very small, the single-lepton + ft channel will show his true potential 
only after the LHC has collected enough statistics. We find that with 300 fb -1 (3000 fb -1 ) 
of integrated luminosity, fix + it searches should be able to exclude scales A s t p > 190 GeV 
(A s,p 210 GeV) if DM is light. Notice that the ATLAS sensitivity study [50] observed a 
quite similar luminosity dependence of mass limits in the case of top-squark searches. 
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Figure 4. The quantities As (left panel) and A p (right panel) as a function of the DM mass. 
The solid lines indicate the results that apply in the case of the 8 TeV mono-jet searches, while the 
dashed curves correspond to our future projections based on 25 fb -1 of 14 TeV data. 


3.4 Infinite top-quark mass limit 

The mono-jet limits on A s,p have been derived in the last two sections by employing the 
exact results for the one-loop pp —> + j amplitudes, but integrating out the mediators 

that induce the DM-SM interactions. In the following, we compare these bounds to those 
that one obtains in the limit of infinite top-quark mass. We call the latter limits (A s,p) m 
and define 

(A5 p) 

Ag p = - S V m ‘- >0 ° -1, (3.13) 

A s,p 

which is a measure of the (in)accuracy of the heavy top-quark approximation. Note that the 
quantities A s,p depend sensitively on the and pr,j cuts imposed in the experimental 
analysis. The two panels in Figure 4 show our results for the quantities A 5 (left) and 
Ap (right). As indicated by the solid curves, for the existing mono-jet searches we find 
A s £ [25,80]% (Ap E [30,90]%). Recalling that the ratio of the $ T + j cross sections 
scales as (l + As,p) 6 , it follows that the nit —> 00 limit overestimates the exact results 
by a factor 4 (5) for small DM mass and that the quality of the approximation rapidly 
degrades with m x , resulting in errors of up to a factor of 32 (48) for the operator Og 
(Op). This clearly shows that in order to infer faithful bounds on the DM top-quark 
contact operators (2.5), one has to calculate the mono-jet cross section keeping the full 
top-quark mass dependence [ 6 , 9-12], Notice that at the 14 TeV LHC the mt 00 limit 
is an even worse approximation than in the LHC run-1 environment. Numerically, we 
obtain A 5 £ [65,105]% (Ap £ [70,120]%), which translates into factors of 19 to 72 (27 
to 108) at the level of cross sections. The observed differences between the 8 TeV and 
the 14 TeV results are easy to understand by remembering that at higher energies the 
imposed and pnj cuts have to be harsher (see (3.1) and (3.9)) to differentiate signal 
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from background. High-energetic initial-state and/or final-state particles are however able 
to resolve the structure of the top-quark loops that generate the 4 - j signal, so that 
removing the top quark as an active degree of freedom becomes less and less justified the 
more stringent the and Ptj selection requirements are. 

The infinite top-quark mass limit is also expected to fail badly at the NLO level. This 
means that taking a K factor obtained from Higgs plus jet production to upgrade the 
LOPS mono-jet cross sections to the NLO level (as done in [11]) is hard to defend from a 
theoretical point of view. A further complication in estimating the size of NLO contributions 
t° o'fid (pp —;► + j ) arises from the fact that in the LHC mono-jet analyses a jet veto is 

imposed. Such a jet veto tends to decrease the importance of fixed-order NLO corrections. 
For instance, in the case of the operator Oq it was found in [22] that the K factor is 
reduced from K ~ 1.5 at fixed order to K ~ 1.1 after including PS and hadronisation 
effects. Although one naively would expect to find a reduction of similar size also in the 
case considered here, we believe that in order to make a definite statement about the size 
of NLO corrections to the $ T + j signal associated to the operators Og p , an exact 0(af) 
calculation of the top-quark loop-induced mono-jet cross section is unavoidable. 

3.5 Comparison of present DM constraints 

Below we present the current LHC exclusion limits on the parameters entering the simpli¬ 
fied model (2.1), comparing them to those stemming from the first LUX results [51] on the 
SI cross section, the latest Fermi-LAT bounds [52] on the velocity-averaged total DM anni¬ 
hilation cross section and the requirement not to overdose the Universe. Our exact collider 
bounds will furthermore be contrasted with the constraints that derive by employing the 
EFT framework (2.5). 

Limits from mono-jet searches 

The different constraints on the four parameters gg M , rn x aR d A/g that characterise 
our simplified scalar mediator models are summarised in the panels of Figure 5. Here 
and below the widths of the mediators are calculated using (2.2). In physical terms this 
means that we assume that the total decay widths Ts,p are minimal. The shown limits are 
obtained from (3.2) and take into account theoretical uncertainties due to scale variations. 
To be conservative, we include these uncertainties by setting the bounds using the signal 
cross sections calculated for £ = 2 (see Section 3.1). 

Turning our attention to the pl ane (upper left panel), we see that for m x = 

100 GeV and Mg = 300 GeV, the present lower limit (3.2) on the mono-jet cross section 
allows to probe only simplified models with g^ M > 3 and > 0.2 (red contour). This 
finding is in line with the observation made recently in [ 12 ] that current fix + j searches 
are not sensitive to weakly-coupled realisations of (2.1). Although a direct comparison is 
difficult, our conclusions also seem to agree with the results presented in [11]. We have 
explicitly verified that even for light scalar mediators with Mg < 100 GeV and m x < 
50 GeV, values of 5^ dm < 1 remain inaccessible, if the mediator width Tg is calculated 
using the formulas (2.2). Our finding that existing mono-jet searches are only sensitive 
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Figure 5. Present mono-jet exclusion regions at 95% CL (red contours) for scalar mediators. In 
the gf M -<7D M plane (upper left panel) the values m x = 100 GeV and Ms = 300 GeV have been 
employed, while in the gg-m x plane (upper right panel) we have identified gs = <?f M = 9dm an( l 
fixed the scalar mediator mass to Mg = 300 GeV. The results in the gs -Mg plane (lower left panel) 
use the same identification and a DM mass of m x = 100 GeV, whereas in the Mg-m x plane (lower 
right panel) the couplings have been set to gg M = gj=) M = 4. For comparison the regions with 
Tg > Mg (brown contours), the current LUX 90% CL constraint on cr^ (solid blue curves), the 
parameter spaces with f l x h 2 < 0.11 (dot-dashed purple curves), the EFT limits (dashed red curves) 
and the regions with Mg > 2 m x (dotted black lines) have been indicated. 


to scenarios with couplings g^ M of order of a few is hence robust against variations of the 
remaining parameters. 

The exact 95% CL exclusion region should be contrasted with the limits that follow 
from an EFT interpretation (dashed red curve) of the + j searches. We observe that for 
the considered values of m x and Mg, the EFT bounds are too strong (weak) for > 2 


- 15 - 
































m x = 100 GeV, Mp = 300 GeV 




9p — 9sm — 9dm i Mp — 300 GeV 



P Pa 

9 sm = 3dm = 4 



Figure 6. Present mono-jet exclusion regions at 95% CL for pseudo-scalar mediators. The latest 
Fermi-LAT 95% CL bound on the total velocity-averaged DM annihilation cross section (<jv x ) is 
indicated by the solid blue curves. Apart from this the same colour coding and choice of parameters 
as in Figure 5 is adopted. 


(•9dm 1$ 2). This feature is easy to understand by noticing that for our choice of parameters 
one has Ms < 2 mt, which implies that Tg ~ T(S I — > xx) c* (.9 dm) The width of the 
scalar mediator thus grows quadratically with and for > 6 one ends up in the 
unphysical situation where Ts > Ms (brown contour). For a broad (narrow) resonance it is 
however known (see for instance [3, 9, 53] for the case of vector and axial-vector mediators) 
that EFT cross sections tend to overestimate (underestimate) the exact results. This is a 
general shortcoming of the EFT framework that can only be overcome by calculating 
signals in a simplified model such as (2.1). 

For comparison we also show in the 3 sm~9dm pl ane the restriction on cjgj provided by 
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LUX (solid blue curve) and the DM relic density (dot-dashed purple curve). Since in the 
case of scalar exchange the elastic DM-nucleon scattering is SI and unsuppressed (cf. (2.8)), 
the limits from the existing direct detection experiments are significantly more stringent 
than the collider bounds, and essentially exclude the entire g^ M DM parameter space for 
Mg = 300 GeV. Notice that the constraints arising from the limits on cr(? T can in principle 
be evaded by assuming that x is n °t stable on cosmological time scales, but lives long 
enough to escape the ATLAS and CMS detectors. We add that the limits from direct 
detection are plagued by systematic errors due to the uncertain local DM density and 
velocity distribution, which play no role in the case of the collider bounds. It is also evident 
that compared to the exact LHC exclusion, the requirement not to have a too high DM 
relic density, i.e. Xl x h 2 < 0.11, further pushes fl’gM DM to larger values. We add that the 
limits following from the relic abundance calculation are more model dependent than the 
remaining bounds, because they depend strongly on the full particle content and all the 
interactions of the underlying theory. For instance, opening up additional DM annihilation 
channels will generically have a more visible impact on Q x h 2 than on o(pp —> + j) and 

Og). These loopholes should be kept clearly in mind when interpreting bounds associated 
to the thermal DM relic density. 

Further insights into the limitations of the EFT description of the mono-jet signal can 
be gained by examining the predictions in the gg-m x (upper right panel), gg-Mg (lower 
left panel) and the Mg-m x planes (lower right panel). A simple criterion that has been 
proposed (cf. [2, 3, 54]) to assess the validity of the EFT approach is to demand that 
Mg > 2 m x . To show how this requirement restricts the domain of the EFT limits we 
have included it into the plots (dotted black lines). From the gs~m x ( gs~Mg ) plane, 
we see that in the region m x > 150 GeV (Mg < 200 GeV) off-shell production of DM 
pairs is numerically important, and as a result the exact exclusion extends far into the 
parameter region with Mg < 2 m x . In the case of the Mg-m x plane, on the other hand, 
one observes that combining the EFT bound with the criterion Mg > 2 m x singles out a 
slice in parameter space that at least qualitatively resembles the real exclusion contour. 
This shows that although simple criteria to gauge the applicability of the EFT cannot be 
used to do precision physics, they are still useful in the sense that they can serve as a sanity 
check of the calculation in the full theory. We finally remark that for the same choice of 
parameterisation of and identical input parameters the shapes and locations of our 
exact mono-jet exclusion contours in the Mg-m x plane do not resemble the results of [12]. 

We now repeat the above exercise for the case of pseudo-scalar s-channel exchange. 
The corresponding results are shown in the panels of Figure 6. One observes that the 
shapes of the exact exclusions (red contours) for pseudo-scalar exchange are quite similar 
to those found in the scalar case, but that the disfavoured regions are larger in all panels. 
The reason for the latter feature is twofold. In the regions Mp > 2 m x (dotted black lines) 
it is a result of (2.7), which implies that for the same model parameters the mono-jet cross 
sections associated to pseudo-scalar exchange are larger than those for scalars by a factor 
of roughly 2. For Mp < 2 m x , one instead has to take into account that the squared 
matrix element of P —>• XX scales as m xx , whereas for the S —> XX channel one obtains 
m xx — 4 m x . The effects of off-shell DM pair production is thus more pronounced if the 
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Figure 7. Exclusion regions at 95% CL for scalar mediators following from the present $ T + it 
searches in the single-lepton channel (green contours). In the 5sm~5dm plane (upper left panel) the 
values m x = 100 GeV and Ms = 300 GeV have been used, while in the gs~m x plane (upper right 
panel) we have set gs = g§ M = 5j=) M and Mg = 300 GeV. The results in the gg-Mg plane (lower 
left panel) use the same couplings and m x = 100 GeV, whereas in the Mg-m x plane (lower right 
panel) the couplings have been fixed to = 4. The regions with Tg > Mg (brown 

contours), the current LUX constraint on cr^ (solid blue curves), the parameter spaces with f l x h 2 < 
0.11 (dot-dashed purple curves), the EFT limits (dashed green curves) and the regions with Mg > 
2 m x (dotted black lines) are also shown. 


mediator couples to X75X rather than to XX- Chir results for the mono-jet exclusion regions 
are again in qualitative (though not quantitative) agreement with [12]. Two other visible 
differences are that the regions (brown contours) where the mediator width exceeds its 
mass, i.e. Tp > Mp , cover more of the shown parameter space, and that direct detection 
does not provide relevant constraints, because the DM-nucleon scattering cross section is 
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Figure 8. Exclusion regions at 95% CL for pseudo-scalar mediators following from the current 
$ T + it searches in the single-lepton final state. The present Fermi-LAT 95% CL limit on the total 
velocity-averaged DM annihilation cross section ( av x ) is indicated by the solid blue curves. The 
colour coding and choice of parameters otherwise resembles the one used in Figure 7. 


spin-dependent and momentum-suppressed for pseudo-scalar interactions. The leading non¬ 
collider constraints arise therefore from the Fermi-LAT measurements of the 7 -ray flux of 
dwarf spheroidal satellite galaxies of the Milky Way (solid blue curves). In fact, one observes 
that combining the LHC constraints with those stemming from ( av x }, Q x h 2 and Tp > Mp 
restricts the allowed parameter space visibly (white regions between the solid blue and dot- 
dashed purple curves). Notice however that the bounds that follow from determinations 
of the velocity-averaged total DM annihilation cross sections are model dependent. For 
instance, they can be weakened significantly in the region m x < mt, if the pseudo-scalar 
mediator does not couple to down-type quarks, i.e. by choosing = 0 and 7 ^ 0 (see 
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remark below (2.1)). Such a choice is compatible with the MFV hypothesis and will not 
affect the mono-jet limits that result from top-quark loops. 

Limits from single-lepton + it channel 

The restrictions on the simplified scalar and pseudo-scalar mediator models that follow from 
the bound (3.6) on the fiducial cross section of pp —> $ T +jbl are shown in Figure 7 and 8, re¬ 
spectively. One first observes that the constraints on the four-dimensional parameter space 
are in the case of the fix -f- tt signal very similar for scalar and pseudo-scalar interactions. 
This is a consequence of the observation made earlier that a (pp — > S — > fix +tt) — a (pp —> 
P fi T + tt) if Ms = Mp, g^ M = gg M , g§ M = g^ M and m x is sufficiently light. From the 
different panels one can furthermore see that for scalar interactions the single-lepton fix + tt 
channel can provide stronger constraints than the fix + j searches in regions of parameter 
space with dominant off-shell production, while for pseudo-scalar mediators this is not the 
case. Given the discussion of the EFT limits in Section 3.2, this is an unexpected finding. 
Qualitatively our observation can be explained as follows. First, in the case of the mono-jet 
signal there is a strong form-factor suppression at work that originates from the momentum 
dependence of the top-quark loop amplitudes. As a consequence, the EFT limits typically 
tend to be too strong when compared to the exact exclusions (see e.g. the Ms~m x plane in 
Figure 5). In the case of the single-lepton fix + tt channel, on the other hand, one can also 
find parameter regions where the opposite behaviour is observed (see e.g. Ms~m x plane in 
Figure 7). This example shows clearly that depending on the dynamic of the considered fix 
process an EFT description can lead to both too aggressive and too conservative bounds. In 
order to determine the exact exclusions a calculation in a simplified DM model is therefore 
mandatory. 

Another interesting characteristic of Figures 5 to 8 is that even so the areas of the 
excluded regions differ their shapes are quite similar, if one considers the same parameter 
plane. This feature can be understood by realising that in their present form the fix + j 
and fix + tt searches are simple cut-and-count experiments that measure the total number 
of events in the tails of distributions such as the fix spectrum. Yet, the size of these tails 
depends to first approximation only on the overall production rate, but is rather insensitive 
to the precise form of the DM-SM interactions that lead to a given final state. This implies 
that while the existing fix searches are well suited to bound/discover DM, they are unlike 
to provide enough information to determine further DM properties. For instance, with the 
existing fix + tt searches it is impossible to distinguish a fix signal arising from S XX or 
P -X- XX- Some of these limitations can however be overcome by studying two-particle (or 
multi-particle) correlations in processes involving fix [10, 55, 56]. In the case of fix + tt 
production, a sensitive probe of the Lorentz structure of the DM-SM interactions is provided 
by the pseudo-rapidity difference A rjb^ (A rji+i~) between the two 6-jets (charged leptons) 
that result from the top-quark decays [18]. Like in the case of Higgs physics (see e.g. [57, 58]) 
studies of the correlations of the top-quark decay products in associated production also 
offer in the context of a fix + it signal unique opportunities to probe the DM mediator 
top-quark interactions. Any dedicated effort at LHC run-2 in this direction is thus more 
than welcome. 
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Figure 9. Exclusion contours at 95% CL for scalar mediators following from hypothetical mea¬ 
surements of a + j signal (red regions) and studies of the $ T + tt single-lepton channel (green 
regions). The shown predictions assume 25 fb -1 of 14TeV LHC data. In the <7 sm~5dm plane (upper 
left panel), we employed rn x = 100 GeV and Mg = 300 GeV, while in the gs~m x plane (upper right 
panel) we have set gg = = g^ M and Mg = 300 GeV. The results in the gg-Mg plane (lower 

left panel) correspond to the same couplings and m x = 100 GeV, while in the Mg-m x plane (lower 
right panel) the couplings have been fixed to g$ M = gj=) M = 3. The regions with Tg > Mg (brown 
contours) and the regions with Mg > 2 m x (dotted black lines) have been indicated for comparison. 


3.6 Future sensitivities of DM searches 

We finally study how the bounds on the parameter space of the simplified scalar and pseudo¬ 
scalar models may improve at future LHC runs. As a baseline for our analysis, we con¬ 
sider 25 fb” 1 of integrated luminosity collected at 14TeV, which corresponds to around one 
year of data taking. The panels in Figure 9 and 10 show our results for s-channel scalar and 


- 21 














9p — 9sm ~ 9dm i Mp — 300 GeV 



9p — 9s m — 9dm 5 m x ~ 100 GeV 



p — p — Q 
9 sm — 5dm — 6 



M P [GeV] 


Figure 10. Exclusion contours at 95% CL for pseudo-scalar mediators following from hypothetical 
measurements of a f T +j signal (red regions) and studies of the $ T +tt single-lepton channel (green 
regions). The shown predictions assume 25 ftC 1 of 14TeV LHC data. The colour coding and choice 
of parameters is identical to the one used in Figure 9. 


pseudo-scalar mediators, respectively. In all plots the red and green contours correspond 
to the 95% CL exclusions obtained from (3.10) in the case of the fix + j signal and (3.12) 
for what concerns the $x + tt single-lepton signature. To allow for an easy comparison be¬ 
tween the present and the future constraints, we have employed in the 9sm~9dmi 9s,P~ m x 
and the gs,P~Ms,P planes the same choice of parameters that has been previously used in 
Figures 5 to 8. One observes that while the shapes of the contours remain qualitatively the 
same, the bounds that one might be able to set with upcoming data will improve notable 
compared to the limits obtained at 8TeV. One however also sees that even with 25fb~ 1 
of 14 TeV data, only model realisations in which the mediators have masses not too far above 
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the weak scale, i.e. Mg p <C 1 TeV, and couple strong enough to the SM, i.e. g s ^ > 1, can 
be explored. This feature is further illustrated by the exclusion contours in the Ms~m x and 
Mp-m x planes, which have been obtained for g§ M = gj§i M = 3 and g^ M = = 3. Notice 

finally that in our sensitivity study the constraints from the single-lepton $ T + it chan¬ 
nel are in the entire parameter and theory space weaker than the restrictions that derive 
from the + j signal. As already mentioned at the end of Section 3.3, since for realistic 
cuts (Jfid(,PP —> $t + tt (—> jblu )) is much smaller than a^d^PP ~> $t + j) ; the $p + tt 
channel will only become competitive to the mono-jet signature at the phase-1 and phase-2 
upgrades at 14 TeV. Realising that the existing $p + tt (—> jblv ) analyses are all recasts of 
top-squark searches (e.g. [16] relies on [43]), the LHC reach might even be improved further 
by trying to optimise these searches to the specific topology of the $p + it signature arising 
in simplified scalar and pseudo-scalar models. This issue deserves additional studies. 

4 Conclusions 

Dedicated searches for DM candidates represent an integral part of the physics programme 
at the LHC. Given our ignorance of the dynamics that may connect the SM to the dark 
sector, it is important that these searches are as model independent as possible and sensitive 
to many different types of DM pair production. One way to achieve this is to employ an 
EFT in which the SM couples to DM via contact interactions. In fact, the EFT approach 
has proven to be useful in the analysis of LHC data, because it allows to derive stringent 
bounds on effective DM-SM interactions involving light and heavy quarks, gluons, photons, 
electroweak gauge bosons and the Higgs that can be easily compared to the limits of direct 
and indirect DM searches. 

In this article, we have studied ftp signatures that result from interactions between 
DM and top quarks. Scalar and pseudo-scalar couplings of this type have been searched 
for in two ways at the LHC. First, by looking for a pp —> $p + j signal, where the DM 
pair is emitted from a top-quark loop, and second by trying to detect the top-quark decay 
products that arise from the tree-level transition pp —>• fip+tt. Our EFT analysis shows that 
the strongest constraints that the 8 TeV LHC can place on both scalar and pseudo-scalar 
DM top-quark interactions come from mono-jet searches. Let us add that DM top-quark 
couplings can also be probed at the LHC in associated production of a single top quark and 
a W boson and t-channel single-top production. We plan to return to these complementary 
search strategies in a forthcoming publication. 

Already at 8 TeV the LHC however operates in a regime where the EFT interpretations 
of the $p signals often do not apply, and these limitations will become more severe at 
a future energy upgrade to 14 TeV. In order to derive faithful bounds on the DM-SM 
interactions, one has thus turn to simplified models, where the dynamics of the contact 
terms is resolved. Motivated by this general observation, we have investigated in detail the 
phenomenology of the simplified models that give rise to interactions between DM and top- 
quark pairs through the s-channel exchange of a scalar or pseudo-scalar resonance. Under 
the assumption of MFV, the parameter space of these simplified models is four dimensional 
and consists of the two couplings (<7 dm) and (<?g as well as the DM and the 
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mediator masses m x and Ms ( Mp ). In our analysis the mediator width T^ (Tp) is not 
treated as a free parameter, but calculated from the remaining parameters. By scanning 
the four-dimensional parameter spaces, we observe that while the signal cross sections 
dependent strongly on f/ D ! M gM , m x and Mg,p . the ratio between the LOPS and the LO 
cross section is, for a given signal region, essentially independent of the choice of model 
parameters. To give an example, in the case of the 8TeV (14TeV) mono-jet search, we 
find a flat ratio close to 60% (45%) (see Appendix A for more details). While we have not 
made use of this feature in our work, this property can be exploited to efficiently generate 
large MC samples with our POWHEG BOX implementations. 

One main outcome of our comprehensive analysis of the simplified models is that even 
with the full LHC run-1 data set theories with couplings g below 1 cannot be tested. By 
increasing the LHC centre-of-mass energy to 14TeV, larger parts of the parameter spaces 
can be explored, but discoveries are still only possible if the mediators have masses of the 
order of the weak scale and couple sufficiently strong to top-quark pairs. Particles with 
such properties contribute indirectly also to other observables such as the total it cross 
section, the Peskin-Takeuchi parameters and necessarily change the properties of the Higgs 
boson, which provides further avenues to search for them. We find that while mono-jet 
searches typically provide the dominant restrictions, in the case of simplified models with 
sizeable couplings between the scalar mediator and top quarks, top-quark pair production 
in association with can allow to better probe parameter regions with dominant off-shell 
production. In contrast, fix + j searches are superior to ftp + tt searches for pseudo-scalar 
s-channel exchange in the entire parameter space. This observation shows the complemen¬ 
tarity of these two different channels and underpins the importance of performing dedicated 
searches for pp —> $p+tt (and likewise pp —> ]fcp+bb) in the upcoming LHC runs. Our exact 
calculations furthermore allow us to determine under which circumstances an EFT interpre¬ 
tation of LHC bounds is possible. While it turns out that simple criteria like Ms t p > 2 m x 
combined with the EFT bounds typically fail to reproduce the exact exclusions, they still 
provide enough information to get a first idea about the sensitivity of different $ T chan¬ 
nels, which makes the EFT framework a particularly useful tool when designing new search 
strategies. We finally compared our exact LHC results to the bounds obtained from direct 
and indirect detection experiments as well as the constraints arising from the thermal DM 
relic abundance and discussed some of the caveats of such a comparison. 

With the start of LHC run-2, collider searches for $p signatures are soon to explore new 
territory, and the large statistics expected at the phase-1 and phase-2 upgrades at 14TeV 
have the potential to revolutionise our understanding of DM. New theoretical developments 
that allow for a better description of both signals and backgrounds have to go along with 
the experimental advances in order to exploit the full physics potential of the LHC. Studies 
based on the simplified DM models we have discussed here may play a key role in this effort. 

Note added 

After communications with us, the authors of [12] revised their mono-jet analysis. They 
found that in their original study they forgot to implement the jet veto, which is supposed 
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to reject events if they contain more than two jets. In the scalar and pseudo-scalar case 
this omission results in fiducial cross sections that are too large by at least a factor of 3. 
After correcting this mistake the results [59] now seem to be in fair agreement with our 
findings. Still at close inspections quantitative differences can be observed. For instance, 
the fact that at 8TeV the exclusion contours in the Ms p—m x planes as shown in [59] do 
not extend up to Mg p = 2 m x calls in our opinion for an explanation. 
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A Generation of mono-jet samples on a budget 

The studies of the simplified DM models performed in this work require to scan the four¬ 
dimensional parameter spaces spanned by , g^, m x and Mg,p■ If the parameter 
space is sampled brute force, i.e. by calculating the LOPS fiducial cross sections for each 
parameter point, this is a time-consuming task that can be done in a finite amount of 
time only on a computer farm. It is therefore worthwhile to ask if accurate results for the 
fiducial cross sections can be obtained without running a full MC chain including a PS 
and a detector simulation. In fact, it is possible to achieve precise results in an efficient 
way, if one makes use of the observation that the ratio of the LOPS and the LO results for 
the fiducial cross sections (which is a measure of the acceptance) of simple $p signatures 
are to very good approximation independent of the specific choice of model parameters. 
This feature is illustrated in the two panels of Figure 11 for the case of a mono-jet signal 
arising from our simplified scalar and pseudo-scalar models. On the left-hand (right-hand) 
side we show the ratio of the LO and the LOPS cross sections in the Ms~m x ( Mp—m x ) 
for ftp + j production at the 8TeV (14TeV). The signal regions are specified in (3.1) 
and (3.9), respectively. Both panels demonstrate clearly that the ratio between the LOPS 
and the LO cross sections is essentially flat in parameter space and only depends on the 
experimental environment, i.e. the centre-of-mass energy of the collisions and the event 
selection requirements that define the fiducial signal regions. Numerically, we find that 
for the considered searches and the imposed generation cuts the cross section ratios lie in 
the range [0.55,0.62] and [0.40,0.46], which corresponds to a relative error of ±6% and 
±7%. Compared to the scale ambiguities that amount to around ±40% in the case of 
mono-jet searches the latter uncertainties are hence subleading and can be neglected to 
first approximation. We have also verified that the above observation applies to the case of 
a mono-jet signature resulting from vector and axial-vector mediator s-channel exchange, 
which can be simulated at NLO and NLOPS level using the P0WHEG BOX implementation 
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Figure 11. Contour plot of the ratio of the LOPS and LO cross sections for mono-jet production 
as a function of the mediator mass Mg ( Mp) and the DM mass m x . The left (right) panel shows 
the results for a scalar (pseudo-scalar) mediator at 8TeV (14TeV). 


presented in [22], Notice that these findings strongly suggest that in the mono-jet case, the 
signal acceptance is a rather flat function in both parameter and theory space. 

From the above discussion it should be clear that to achieve accurate results it is 
sufficient to calculate the ratio between the fiducial cross sections before and after including 
a PS and a detector simulation for a few parameter points only. The determined ratio can 
then be used to promote the fixed-order results to the true fiducial cross sections. If this 
is done the problem of calculating the LOPS (or NLOPS) results boils down to generating 
the LO (or NLO) fixed-order fiducial cross sections efficiently. In the case of the mono¬ 
jet signal this can be done by generating Born-level configurations with a suitable cut on 
the minimal PT.j 1 of the leading jet. Rather than imposing a fixed cut and generating 
unweighted events it turns out to be advantages [22, 60] to generate weighted events of 
relative weight w = {Pt xx + p't sup)/pT xx - Here PT,xx H ie transverse momentum of the 
DM pair and pr sup is a parameter that suppresses the generation of phase space points with 
small Pt,xx • Hi th e case °f rno no-jet searches one has to choose px ,sup sufficiently below 
the fix restriction imposed in the experimental analysis. Using this method to populate 
the phase space will lead to a rather uniform distribution of events in the entire Px, X \ range, 
but the LOPS (or NLOPS) accuracy is recovered, because the few events at low Px, xx 
have a large relative weight w. In this way, the correct cross section is reproduced. To allow 
for maximal flexibility both methods of sampling the phase space have been implemented 
in our P0WHEG BOX add-on that can simulate the signal arising from top-quark loops. 
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